x=[0.0:4.0:0.2]
y=sin(5*x)
ey=replicate(0.1,size(x,/n_elements))
ex=ey
a=errorplot(x,y,ex,ey)
end